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within the critical state approximation. The time dependences of the current density 
profiles, magnetic flux lines, local power dissipation and magnetic moment are fea- 
tured. Noticeable non-homogeneous dissipation and field distortions are displayed. 
Also, significant differences between the obtained AC-losses and those predicted by 
regular approximation formulas are reported. Finally, an outstanding low pass filter- 
ing effect intrinsic to the magnetic response of the system is described. 
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The practical configurations of type-II superconducting wires exhibit complicated non- 
linear and hysteretic behavior under oscillating electromagnetic fields. Thus, related to the 
design of power applications with these materials, considerable efforts have been made to 
comprehend the factors determining the AC losses under typical operating conditions. Ma- 
jor features of the macroscopic electromagnetic behavior have been captured by Bean in the 
phenomenological so-called critical state model (CSM).— According to this, magnetization 
currents are induced at the periphery of the superconductor when external flux variations 
occur. Such currents distribute across the section of the sample with a density equal to 
so-called critical value at a given temperature and field, J c . Although simple in idealized 
configurations, the Maxwell equation problem arising from the CSM statement becomes 
awkward when realistic systems are considered. 

In this letter, we consider the coupling between simultaneous oscillating sources acting on 
a superconducting wire. In particular, the simultaneous action of an AC transport current 
and a transverse magnetic field will be studied. Notice that such a configuration is a basic 
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FIG. 1. (Color online) Sketch of one of the analized cases. Here, B™ ax = 2, Jf r ealt = 0.75J C , and 
the time-step corresponds to B y = —0.64, I tr = —0.24. Units are (hq/A-k)J c R for B, irR 2 J c for I tr , 
J C R ?J for M, and (/xo/4tt) J^R 4 for E-J. The experimental process (B y ,It r ) is summarized in the 
central plot, and the physics behind this kind of system can be easily understood analyzing the 
lattachcd videol 
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model for coil systems in which each wire is under the action of its neighbors. Apparently, the 
problem is 2D in nature and the main complication is as follows. When the external sources 
change in time, a sort of free-boundary problem has to be solved. Thus, the magnetic flux 
density variations penetrate from the sample's surface and an evolutionary flux-front profile 
is defined. Determining such front (or the core within) is a main mathematical challenge, 
and has been tackled by several methods.- - — As an important handicap, in many cases as 
the one studied here, the core does not remain static along the electromagnetic time- varying 
process, even for the case of a fully penetrated sample. It is apparent that the widely used 
concept of an effective radius description for the free-boundary can fail markedly if the core 
shifts within the sample (as may be observed in FigJT]) and a center of mass is not well 
defined. As a consequence the so-called front tracking methods cannot easily guarantee a 
unique physical solution for the problem. 

Here, we adopt the most popular trend in the analysis of electromagnetic applications, 
i.e.: numerical simulations implementing finite-element methods. As an advantage, such 
techniques work without explicit inclusion of the free boundary. The whole superconducting 
region is involved in the calculation and the boundary is obtained as a part of the solution. 
Taking advantage of this methodology, we have performed a systematic investigation of the 
electromagnetic response of the superconducting wire. Some outstanding predictions have 
been obtained and are reported here. In particular, (i) we feature the relevance of localized 
power dissipation within the sample, (ii) as regards the averaged physical quantities, we 
will show that some standard approximations for the power dissipation per cycle have to be 
revised, and finally (iii) an intriguing low pass filtering effect is announced. 

Going into detail about our calculations, we use a discrete formulation that solves Fara- 
day's law iteratively in a mesh of circuits that carry the macroscopic electric current.- Cal- 
culations are performed under the material law restriction for the current density, that in 
this case reads |J| < J c . Equivalently, one may use conductivity law E = pj with p(J) = 
when | J | < J c and p(J) — > oo if |J| > J c . 

Under the assumption of translational symmetry, the finite-element implementation no- 
ticeably simplifies. Explicitly, we discretize the sample's cross-section by a collection of 
straight infinite elementary wires, each of them carrying a current Ii = JiSi with Jj the cur- 
rent density and Si the cross sectional area of the element. As it has been shown in previous 
work,- a variational formulation that suits the discrete modeling is possible for our electro- 
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magnetic problem. Thus, one can show that in quasi-steady regime (excellent approximation 
for the large scale application frequencies) the discrete form of Faraday's law is obtained by 
minimizing a magneto- quasi- static field Lagrangian (of density C = (B* +5 * — B*) 2 /2), cou- 
pling successive time layers and under prescribed sources and material law. In our case, and 
by using standard electromagnetic manipulations, the quantity to be minimized becomes 



with Ij the set of unknown currents for the collection of elements, My the mutual induc- 



A e the vector potential related to external sources (applied transverse field). Optimization 
has to be performed under the restriction of applied transport current in the cross section Q, 
Sign = I*n an d f° r the critical state law Li < I c . This has been done by using specialized 
large scale constrained minimization algorithms as discussed before.- 

Some technical comments are worth of mention for the case under study. In general, the 
electromagnetic manipulations leading to Eq.fll]) may be done as follows. By using that one 
can split the cross section of the cylinder into a high number of elements (wires), it is justified 
to assume that the elementary inner potential is created by a uniform current density and 
contributes to the energy as a constant, thus one can use Ma = Hq/A in Eq.flJJ). On the 
other hand, using the logarithmic expression for the two-dimensional Green's function for 
the outer potential of the wires one gets M™ 1 = (/io/47r)ln(rjj) for the mutual inductances. 
These expressions have been obtained under the condition of continuity for the magnetic 
potential but are arbitrary save to a global constant (gauge invariance). In the absence of 
transport current, one can show that such constant may be obviated in the theory. The 
reason is that, when deriving Eq.flT]), spatially constant terms from the potentials of the 
wires are multiplied by X^ii (= 0). However, for problems with transport, unless one 
cares about such terms, some calculations, as the value of E may be tampered. In order 
to show how this arises, we recall that, generally speaking, physically acceptable electric 
fields have to be expressed in the form E = — d t A — V0, including an electrostatic like term. 
For long wires, it can be argued that V0 has to be constant in space, and thus one has 
E = — [d t A — C(t)]k = [<9iv4']k where A' works as a calibrated potential.— Taking advantage 
of the fact that arbitrariness is only up to a constant, the situation may be tackled by 
progressively determining C{t) according to the physical condition E = at those points 
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(1) 



tance matrix between such elements i and j, Ji 



the solution at the previous time layer, and 



where the magnetic flux does not vary. 

Simulations have been performed for the triangular oscillating process displayed in Fig. [TJ 
The following quantities have been focused: (i) magnetic field lines derived from B = VxA, 
(ii) the sample's magnetic moment (per unit length) M = (1/2/) f n r x J, (iii) the local 
density of power dissipation (E • J), and the hysteretic AC losses per unit time and volume 
for cyclic excitations of frequency u calculated as L = (u/2ti 2 R 2 ) ^y c dt j n E ■ J. Here, fx. 
denotes a full cycle of the time- varying electromagnetic sources. Our results, and conclusions 
are developed along the following paragraphs. 

a. Flux penetration profiles. Fig. [1] shows some of the results obtained for one of the 
experimental processes considered along this work. The dynamics of the electromagnetic 
quantities can be followed in detail by means of the attached video. Several aspects to be 
noticed are: (i) the position of the flux free region (core) and the current density profiles are 
not axially symmetric by consumption of the magnetization currents. Related to this, (ii) 
strong distortions of the magnetic flux lines appear around the superconducting cable. This 
is especially notorious when the value of the applied field and the transport current approach 
zero, (iii) Outstandingly, the local density of power dissipated along the process also displays 
a strong localization. In fact, a higher heat production (and transfer for constant T) is always 
predicted for half cross-section of the superconducting wire. We argue that this phenomenon 
could increase the quench probability. Furthermore, the proper determination of the active 
area depends on the history of the first branch of the experimental process, e.g. in Fig. [1] 
a positive slope in both B y (t) and hrif) determines maximal power dissipation towards the 
positive x-axis. 

b. Analysis of AC losses. Fig. [2] shows the calculated variation of the hysteretic AC 
losses for long wires in terms of the amplitude of the coupled electromagnetic sources, 
L(B™ ax , It^ ak )- Our results are compared to those obtained from several analytical ap- 
proximations^ customarily applied for non-coupled periodic sources (see caption at Fig. [2]) 
and Bean-like superconductors (J c = constant). The main significance of our results is 
that linear superposition only makes sense for high magnetic fields and moderate (or low) 
currents. This complements previous work on the rectangular geometry^ and adds new 
perspectives on the validity of approximation formulas. We emphasize the failure of assum- 
ing simple linear superposition of the classical formulas.- In fact, linear approximations as 
L(B) + L(Itr^) and L(B, ( see Ref- [^) can notoriously underestimate or overestimate 
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FIG. 2. (Color online) Calculated hysteretic ac losses per cycle for simultaneous oscillating trans- 
port current and magnetic flux density, L(B,I^ C ^). The different AC transport conditions are 
denned by the corresponding field amplitudes and labelled peak currents. Our results are directly 
compared with the conventional approaches^ for: (Left) isolated sources L(B) and L{I^ c} ), and 
the intuitive linear superposition of both of them. (Middle) oscillating magnetic flux density under 
a constant transport current condition L(B, (Right) The whole set of results are also plotted 

in linear scale. Losses units are are 



the real losses. Notice that, even for the best situation, L(B) + L(// r ), at low amplitudes 
of B the differences between predicted losses can overcome 100% (e.g. at L(l, < 0.25)). 
Then, they decreases as the amplitude of ij;^ increases (~ 15% for 4 ac} ~ 1). On the 
other hand, if the amplitude of B is as high as B p = 8 (the magnetic full penetration field 
for zero transport current) the difference can oscillate between ~ 20% (l{ r ac} ~ 0.25) and 
~ 55% (Itr^ ~ 1) or even higher for ~ 0.75. 

On average for the extensive set of experiments analyzed, the highest differences (~ 
100%) are observed around the condition (B p /2, I c /2). Specially at low values of B, the 
differences between analytical approaches and our numerical calculation can overcome 100%. 
In conclusion, for a proper determination of the hysteretic AC losses in systems with coupled 
electromagnetic sources (£>, // r ac ^) somehow sophisticated analysis resources are needed, even 
for relatively simple configurations as the one studied here. 

c. Magnetic moment cycles (low pass filtering) Fig. [3] shows the dynamics of the mag- 
netic moment component along the external field M y in terms of the evolution of the mag- 
netic sources B y (t), I tr (t). Notice that, only for small values of the applied transport current, 
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FIG. 3. (Color online) The renormalized magnetic moment M y /Max{M y ) as a function of the 
applied sources, It r (left), B y (middle), and its temporal evolution (right). Several experiments are 
shown according to the maximal intensity of the applied sources as was sketched in Fig. [TJ Indeed, 
each row corresponds to a singular choice of lf^ ak =0.25, 0.5, 0.75, and 1.0. The amplitude of 
the oscillating magnetic field has been assumed also to change to consider the whole spectrum of 
possible experiments. Here we show the corresponding results for B y nax =1, 2, 4, 6, and 8, as can 
be deduced directly from the plots. 

nearly standard Bean-like loops are obtained. In particular, we recall that the flat satura- 
tion behavior of M y for high field values progressively disappears by increasing the value 
If^ ak . Remarkably, this phenomenon ends up with a symmetrization of the loop (both as a 
function of I tr and B y ), that exhibits a characteristic eye-shape. As a consequence of such 
process, an outstanding low pass filtering effect is predicted for the experimental situation 
described in this paper, which may lead to envisage new applications for superconducting 
wires. Thus, a plot of the induced magnetic moment as a function of time (rightmost col- 
umn in Fig. [3]) demonstrates that the triangular input excitation produces a nearly perfect 
sinusoidal output M y (t). 

d. Conclusion. In this work, we have investigated the electromagnetic response of a 
superconducting circular wire under simultaneous AC sources through a numerical imple- 



mentation of Bean's critical state model. The main physical assumptions are an infinitely 
steep E(J) law that goes from zero to infinity at the critical point J c . In simple systems, 
these assumptions are known to lead to hysteretic (rate independent) losses with simple 
relations between monotonic and cyclic quantities. Here, we show that even for the 2D 
configuration studied, one can find intriguing phenomena, unexpected in simplified models 
that are based on plain linear superposition. Thus, we have obtained notorious localization 
effects in the density of power dissipation, strong field distortions, important failures of the 
customary approximation formulas for the bulk quantities, and also predict an outstanding 
low pass filtering effect in the magnetic response. 

Funding of this research is gratefully acknowledged.— 



REFERENCES 

^.P. Bean, Phys. Rev. Lett. 8, 250-253 (1962); Rev. Mod. Phys. 36, 31-39 (1964); 

2 A. Badfa-Majos, C. Lopez, and H. S. Ruiz, Phys. Rev. B 80, 144509 (2009); H. S. Ruiz, 

and A. Badfa-Majos, Supercond. Sci. Technol. 23, 105007 (2010); H. S. Ruiz, C. Lopez, 

and A. Badfa-Majos, Phys. Rev. B 83, 014506 (2011). 
3 M. Ashkin, J. Appl. Phys. 50, 7060 (1979). 
4 M. Ashkin, G.R. Wagner, J. Appl. Phys. 60, 2477 (1986). 

5 W.J. Carr, AC loss and macroscopic theory of superconductors (Gordon & Breech, NY, 
1983). 

6 A.V. Gurevich, R.G. Mints, A.L. Rakhmanov, Physics of composite superconductors, 
Nauka Publishers, 1987 (in Russian) [Engl, translation: Begell House, NY, 1997). 
7 C.Y. Pang, A.M. Campbell, PG. McLaren, IEEE Trans. Mag. 17, 134 (1981). 
8 K.L. Telschow and L.S. Koo, Phys. Rev. B 50, 6932 (1994) 
9 Y.E. Kuzovlev, JETP Lett. 61, 1000 (1995) 

10 F. Gomory, R. Tebano, A. Sanchez, E. Pardo, C. Navau, I. Husek, F. Strycek, and P. 

Kovac, Supercond. Sci. Technol. 15, 1311 (2002) 
n B. ten Haken, J.J. Rabbers, H.H.J, ten Kate, Physica C 377, 156 (2002) 
12 D. Karmakar, K.V. Bhagwat, Physica C 398, 20 (2003) 

13 L. Rostilla, S. Brisigotti, and G. Grasso, J. Supercond. Nov. Magn. 24, 313 (2011). 
14 L. Prigozhin and V. Sokolovsky, Supercond. Sci. Technol. 24, 075012 (2011) 

8 



E. Pardo, F. Gomory, J. Souc, J. M. Ceballos, Supercond. Sci. Technol. 20, 351 (2007) 
Spanish CICyT and FEDER program (project MAT2008-05983-C03-01), the DGA grant 
T12/2011, and Spanish CSIC (JAE program). 



9 



Profiles of Current [ -I (blue) , 1 1 (red) ] Density of power dissipation 




FIG. 4. The video has been uploaded to the Data Conservancy Pilot Project 
of arXiv repositories. Also, this and more videos are available to download in: 
http : //www . unizar . es/departamentos/f isica_mat_condensada/people/hsruizr/ 

Movie Caption: 

Sketch of one of the analized cases along this letter. In this case, B™ ax = 2, and Jf r eak = 0.75J C 
(Figure [I] is a frame of this video). Units are {^q/A-k)J c R for B, ttR 2 J c for It r , J C R? for M, and 
(/jLo/Att) J^R 4 for E-J. Here, we have assumed a quasistationary time-step defined by the experi- 
mental process (B y ,It r ) which is summarized in the central plot. Top-Left. The magnetic field 
lines (projected isolevels of the vector potental over the wire cross-section) and their corresponding 
profiles of current are shown. The consumption of local magnetization currents and the gener- 
ated field distortions around the superconductor wire are directly visualized. From time to time, 
straigth isolines in zones free of current are ploted as a consequence of the number of isolevels which 
has been required to be plotted. They must be understanding only as a visual effect introduced 
by the graphical processing. Top- Right: Dynamics of the density of power dissipation along the 
cross-section of the superconducting wire. A clear non symetric distribution of the heat transfer in 
a cyclic process is always observed. Bottom-Left: The magnetization loop for the superconducting 
wire in terms of the renormalized electromagnetic sources (B y /B™ ax ■ 0.75, 0.75) is shown. Bottom- 
Right: The magnetic moment is shown in terms of the virtual time defined by the dynamics of 
the electromagnetic sources. A low pass filtering effect in the magnetic response as Jt r — > J c is 
predicted. 
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